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The spin-1/2 Ising model on a square lattice, with fluctuating bond interactions 
between nearest neighbors and in the presence of a random magnetic field, is inves- 
tigated within the framework of the effective field theory based on the use of the 

*rt ' differential operator relation. The random field is drawn from the asymmetric and 

anisotropic bimodal probability distribution P(hi) = p5(hi — hi)+q5(hi+chi), where 
the site probabilities p, q take on values within the interval [0, 1] with the constraint 
p + q = 1; hi is the random field variable with strength h\ and c the competition 
parameter, which is the ratio of the strength of the random magnetic field in the two 
principal directions +z and — z; c is considered to be positive resulting in competing 
random fields. The fluctuating bond is drawn from the symmetric but anisotropic 
bimodal probability distribution P(Jij) = \{5{Jij — ( J + A)) + 5(Jij — (J — A))}, 
where J and A represent the average value and standard deviation of Jjj, respec- 

Tj - ! tively. We estimate the transition temperatures, phase diagrams (for various values 

of system's parameters c,p, hi, A), susceptibility, equilibrium equation for magneti- 
zation, which is solved in order to determine the magnetization profile with respect 
to T and h\. 
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1 Introduction 



Prediction of the critical behavior of modified spin models in the presence of 
site or bond dilution, random bonds, random fields has been the subject of 
many studies in the last decades; this modification brings about considerable 
changes in the critical behavior of these systems, such as replacement of a first- 
order phase transition (FOPT) by a second-order phase transition (SOPT), 
depression of tricritical points and critical end points, new critical points and 
universality classes, etc [1,2,3]. The study of the aforementioned disordered 
systems is based on the standard models, such as Ising, Blume-Capel, Baxter- 
Wu, Heisenberg, etc, modified accordingly to meet the requirements under 
consideration. Furthermore, extensions and versions of these models can be 
applied to describe many other different situations, such as multicomponent 
fluids, ternary alloys, 3 He - 4 He mixtures, in addition to the magnetic systems 
for which these were initially conceived. The most extensively investigated 
model in statistical and condensed matter physics is the spin-1/2 Ising model 
(IM), since its two dimensional version, without an external magnetic field, 
was analytically solved by Onsager; it is a simple one relative to other models 
of cooperative phenomena and has a wide range of applicability to real phys- 
ical systems. Its three-dimensional version has not yet been solved exactly, 
for which the only existing results are either from the renormalization group 
calculations or series expansions and are thus considered to be the "exact" 
ones. In its modified versions, it exhibits a variety of multicritical phenomena, 
such as a phase diagram with ordered ferromagnetic and disordered paramag- 
netic phases separated by a transition line that changes from an SOPT to an 
FOPT joined at a tricritical point (TCP); besides these, critical points, critical 
end points, ordered critical points of various orders, re-entrance can appear 
as in the presence of random fields. The multicritical phenomena appear in 
systems presenting competition among distinct types of ordering and there 
are numerous circumstances in which this kind of phenomenon can arise. In 
ferromagnetic systems in the presence of random fields, the competition is be- 
tween the parallel and random ordering, causing, occasionally, the conversion 
of a continuous transition into an FOPT and the subsequent appearance of 
TCP as well as re-entrance in some cases [4]. Random field effects on magnetic 
systems have been systematically studied not only for their own theoretical 
study but for their experimental importance, as well. 

The methods used for the study of the IM are the mean field approximation 
(MFA), high temperature expansions, series expansion, renormalization group, 
Monte Carlo, effective field theory (EFT). The eldest one is the MFA, which 
is very popular because of its simplicity and has played an important role for 
the description of cooperative phenomena for many years; it gives qualitative 
agreement with experimental results for many of the physical quantities in- 
volved in a phase transition. However, MFA has some unsatisfactory results 



in describing correctly the critical region because of the omission of correla- 
tions when its results are compared with those of experiments. MFA is still 
implemented because it is not complicated and can give a general view of 
the expected behavior for the system under consideration. Many efforts were 
attempted towards improving MFA, such as the effective field theory by Hon- 
mura and Kaneyoshi [5]. EFT relies on introducing a differential operator 
into the exact spin correlation function identities obtained by Callen [6] and 
using the van der Waerden spin identities, which improves substantially on 
the standard MFA [7,8]. The EFT approach is general and may be applied 
to systems with any spin value, adapting the van der Waerden identities ac- 
cordingly. This procedure presents a great versatility and has been applied to 
several occasions such as pure, site- and bond-random Ising model, although 
this procedure shall not yield accurate values for the physical quantities in the 
critical region due to the absence of long range fluctuations. EFT has already 
been used in numerous physical problems as a tool to study the magnetic be- 
havior of complex spin systems, such as diluted ferromagnets [9,10,11], pure 
anisotropic systems [12], disordered systems [13,14], cylindrical nanowires [15]. 

The Hamiltonian we shall adopt is that of Ising model with nearest neighbor 
interactions, in which the bond Jy between two neighboring spins varies from 
pair to pair randomly; the system is also under the influence of a random field 
hi, varying from site to site; both random variables are drawn from a suitable 
probability distribution function (PDF). 

The Hamiltonian describing the above system is 

H = - J2 J iA$J ~ E hiSi , Si = ±l (1) 



The summation in the first term extends over all nearest neighbors and is 
denoted by < i,j >; in the second term hi represents a random field that 
couples to the one dimensional spin variable Si. The Hamiltonian describes the 
competition between the long-range order (expressed by the first summation) 
and the random ordering fields. The presence of random fields requires two 
averaging procedures, the usual thermal average, denoted by angular brackets 
(...) and disorder average over the random fields denoted by (...) r for the 
respective PDF's, which are usually a version of the bimodal or Gaussian 
distributions. 

The random exchange integral Jij between the sites % and j is drawn from the 
symmetric and anisotropic bimodal PDF, 

P(Jij) = \Wa ~(J + A)) + KJn -(J- A))} (2) 



where J and A represent the average value and standard deviation of Jy, 
respectively, implying that the exchange integral is fluctuating. This choice of 
the PDF implies that both values (J — A, J + A) of the random exchange 
integral are equally probable. 

The PDF for the random fields hi is also of the bimodal type 

P(h i )=p8(h i -h 1 ) + q8(h i + h 2 ) (3) 

where h 2 = chi, c is the competition parameter and is considered to be positive 
so that the random fields are competitive. The site probability p is the fraction 
of lattice sites having a random magnetic field with strength h±, while the rest 
sites have a field with strength (— h 2 ) with site probability^ such that p+q = 1 
and the usual choice was p — q — |, symmetric case, [4,16,17,18]. As far as 
the PDF (3) is non symmetric, the mean value for hi does not vanish a priori, 
but, instead, is < hi >h= (p — cq)h\, an immediate result is that the system 
is under the influence of a magnetic field in case p ^ q and c^ 1.0. 

One of the main issues was the experimental realization of random fields. Fish- 
man and Aharony [19] showed that the randomly quenched exchange interac- 
tions Ising antiferromagnet in a uniform field H is equivalent to a ferromagnet 
in a random field with the strength of the random field linearly proportional to 
the induced magnetization. This identification gave new impetus to the study 
of the RFIM, the investigation gained further interest and was intensified re- 
sulting in a large number of publications (theoretical, numerical, Monte Carlo 
simulations and experimental) in the last thirty years. Although much effort 
had been invested towards this direction, the only well-established conclusion 
drawn was the existence of a phase transition for d > 3 (d space dimension), 
that is, the critical lower dimension d\ is 2 after a long controversial discussion 
[20], while many other issues are still unanswered; among them is the order 
of the phase transition (first or second order), the universality class and the 
dependence of these points on the form of the random field PDF. The study 
of RFIM has also highlighted another feature of the model, that of tricritical- 
ity and its dependence on the assumed distribution function of the random 
fields. According to the mean field approximation, the choice of the random 
field distribution can lead to a continuous FM/PM boundary as in the single 
Gaussian probability distribution, whereas for the bimodal one this bound- 
ary is divided into two parts, an SOPT branch for high temperatures and an 
FOPT branch for low temperatures separated by a TCP at kT*/(zJ) = 2/3 
and hl/(zJ) = (kT*/(zJ)) x argtanh(l/ v / 3) ~ 0.439 [16], where z is the co- 
ordination number and k the Boltzmann constant, such that for T < T* and 
h > h l c the transition to the FM phase is first order for the symmetric case 
p — \. However, this behavior is not fully elucidated since in the case of the 
three-dimensional RFIM, the high temperature series expansions by Gofman 
et al [21] yielded only continuous transitions for both probability distribu- 



tions, whereas according to Houghton et al [22] both distributions predicted 
the existence of a tricritical point, with h\ = 0.28 ± 0.01 and T c * = 0.49 ± 0.03 
for the bimodal and a* = 0.36 ± 0.01 and T* = 0.36 ± 0.04 for the single 
Gaussian. In the Monte Carlo studies for d = 3, Machta et al [23], using single 
Gaussian distribution, could not reach a definite conclusion concerning the 
nature of the transition, since for some realizations of randomness the mag- 
netization histogram was two-peaked (implying an SOPT) whereas for other 
ones three-peaked, implying an FOPT. 

In this investigation, we study the EFT applied to the random field spin- 1/2 
Ising model on a square lattice with random bond nearest neighbor inter- 
actions, calculating the relevant thermodynamic quantities, such as critical 
temperatures, susceptibility, magnetization profiles with respect to the tem- 
perature T and random field h\. The paper is organized as follows: In the next 
section, the EFT formalism is introduced and applied to the Ising model in 
the presence of a random field and random bonds interactions; the equation 
of state for the magnetization is derived. In section 3, the phase diagrams and 
magnetization profiles are calculated as functions of system's parameters; we 
close with the conclusions in section 4. 

2 Model and formalism 

Considering the Hamiltonian (1), each spin is under the influence of a molec- 
ular field hi 

z 

hi = Y^ JijSj + hi (4) 



where z is the coordination number. For the current model the starting point 
is the Callen exact spin correlation function identity [6] 



{(Si)) r = ((tanh 



P^JijSj + f3hi 



»r (5) 



where the summation takes all the nearest neighboring spin sites of i and << 
• • • >> r indicates the thermal and random configurational averages. Following 
Honmura and Kaneyoshi [5], introducing the differential operator D = -^- into 
(5) and using also the van der Waerden identity 

J^D = cogh ( j.. D ^ + s . sinh ^j.. D ^ ( 6 ) 

Eq. (5) transforms into 



((^))r = ((exp 



3=1 



)) T .(tanh[x + ^/i i ]) x= o 



((J] [cosh(/3 J^D) + Sj sinh(pj tJ D)])) r f(x, /3hi) x=0 (7) 

3=1 



where 

f(x,j3hi) = pt&nh(x + j3hi) + gtanh(x — fih-i) (8) 

Before proceeding to the calculations, we shall consider the following ap- 
proximations to make the problem mathematically tractable, (a) the con- 
figurational average of spins and exchange integral are taken independently, 
{{Sjf(Jij))) — ((Sj))(f(Jij)), (b) the exchange integrals J^ for different j 's 
are also independent of each other, (fi(Jij)f2(Jik)) — (fi(Jij)){f2(Jik))- Under 
these assumptions the relation (7) is written as 

((S t )) r ^J2Q nz J2 ((SnS h ---S jn )) r (9) 

«=1 jl,J2---jn = l 



where all the spins Sj l , Sj 2 , •••, Sj n are nearest neighbors of Si and the coefficient 

^Cnz IS 

Q nz = (coshiPDJ^y^isinh^DJ^rfix, (3h z ) x=0 (10) 



As far as the tanh(rr) is an odd function of its argument, only odd n appear 
in Q nz for the square (z = 4) lattice. Considering the PDF (2) for J^, the 
averages in (10) are 

(cosh(/3 DJij)} = cosh(/3 JD) cosh(/3AD) 

(smh^DJij)) = smh(/3JD) cosh(/3AD) (11) 

thus (10) becomes 

Q nz = cosh z - n ((3JD) sinh n {f3JD) cosh z {f3AD)f{x, Ph t ) x=0 (12) 

which, for even z, takes the form 

z/2 z/2 

Qnz = E h P sinh(2£/3 JD) £ \a ( u z) \ cosh(2u^AD)f(x, ph t ) x=G 
e=i v=o 



z/2 z/2 

E^ n) E l°^l sinh(2£/3J J D) cosh(2vPAD)f(x,Phi) x=0 (13) 

^=1 ^=0 



using the operator relation e wD f(x) = f(x + w), we obtain for the summand 
in (13) 



sinh(2£(3 JD) cosh(2u(3AD)f(x, (3h 



i)x=0 



j{p ( tanh[/3(2£J + 2z/A + fa)] + tanh[/3(2£J - 2z/A + fa)] + 



tanh[/3(2£J - 2z/A - ^)] + tanh[/3(2£J + 2z/A - fc^^ 
g ( tanh[(3(2£J + 2z/A - fa)] tanh[(3(2£J - 2uA - fa)] + 
tanh[/3(2£J - 2z/A + fa)] + tanh[/3(2£J + 2z/A + fc 2 )]) } 



- gb ,(0,J,A,fa,fa) (14) 



thus Eq. (13) is written as 

1 z/2 z/2 

Q™ = 7 E *#° E I4 Z) MA ■/, A, *i, h 2 ) (15) 

4 £=1 i/=0 



with 



a„ 



z/ 



Oi = — 0} 



2 ■ 



bf = u2 f ,. [(2f) 2 -(3,-2)]|af ) | 
lli=ol z "U 

fef = 4 f .s [(2^) 4 - I0(z - 2)(2£) 2 + 15z 2 - 50* + 24]\a < f ) \ (16) 

Applying this procedure to the linear chain (z = 2) without both types of 
randomness, namely, fa — 0, A = 0, (9) implies 

M = Qi 2 ((^i) + (S 2 )) = 2Q 12 (S 1 ) = 2Q 12 M (17) 



where Q\ 2 = |tanh(2/3J) from relations (13), (14), (15), (16), so that for the 
respective critical point we get 2Q\ 2 = 1 or tanh(2/3 c J) = 1 which implies that 

m-j)- 1 = -^ = (is) 



in accordance with the known exact result for the one-dimensional IM, thus 
improving on the MFA result that is (/c.bT c )/J = 2. 

For the disordered square lattice (z = 4), the respective Q nz functions are 



<5i4 — ■^g(2G'i + G2) 

Q U = L(-2G 1 + G 2 ) (19) 



where Gk = 3gko + ^Qki + 9k2 and the g's functions are defined in (14). The 
resulting equation for the equilibrium magnetization from (9) is 

m = AQum + 4Q 34 m 3 (20) 



which admits two solutions, the paramagnetic one, m — 0, and the ferromag- 
netic one given by 




m = ±J -pg (21) 



whereas the critical boundary characterizing the ferromagnetic/paramagnetic 
phases is determined by the condition m = and results as a solution to 
the equation 4Qi4 = 1, which, in case both the random field strength h\ and 
exchange integral deviation A vanish, then it converts into 

2 tanh(2/3 c J) + tanh(4/3 c J) = 2 (22) 



as was also found in [7], and the critical temperature results as kst c /J = 
3.0898 . . ., which is closer to the exact one k B tf act / J = 2/ ln(l + ^2) than the 
MFA one kBt^ FA /J = 4. At this point, it has to be noted that the decoupling 
scheme recalled in Eq. (20) behaves as a better approximation than the one in 
the MFA, since within the EFT framework the kinematics relations are treated 
exactly (a 2 = 1) through the van der Waerden identity; in the present case, 
EFT neglects correlations only between different spin variables, whereas the 
MFA neglects any kind of correlation, namely, the self- and multi-spin ones. 

3 Numerical results. Phase diagrams. Magnetization profiles 

An important feature of the magnetic or fluid systems is their phase diagram, 
the temperature against a suitably chosen parameter, which, in the present 
case, allows the investigation more clearly of the effects of randomness itself 
on the transition temperature of the random system. In the current study, 
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Fig. 1. (Color online) Random field hi variation of the critical temperature T - 
phase diagram (hi — T). In panel (a), for the upper group of graphs c = d = 0.0, 
for the lower group c = 0.0 d = 1.0. In each group from up to down 
p = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0. In panel (b), d is fixed as d = 0.5, c = 0.5 (upper solid 
graphs), and c = 1.5 (lower dashed graphs), labelled by the respective p— values. 
The solid curves are arranged in increasing p from up to down, whereas the dashed 
ones are arranged in increasing p from down to up. In panel (c), c is fixed as c = 0.5, 
d = 0.00001 upper solid graphs, and d = 0.5 lower dashed graphs. The plots in both 
groups are labelled by the respective p— values. T in units of J/k. 

there are various parameters, such as c, A,p, hi, so that we can introduce a 
variety of phase diagrams by plotting the temperature with respect to one 
of these parameters keeping the other ones fixed. However, in the numerical 
calculations to follow the main quantities /3, A, hi, h 2 are written with respect 
to the average value J of the exchange integral, that is, /3 = J/3, d = A/ J, hi = 
hi/ J, h 2 = h 2 /J, scaled quantities; this choice influences Eqs. (2) and (3): the 
PDF in (2) is written as P^) = \{8{Jij - (1 + d)) + 5{J i3 - (1 - d))} so 
that by choosing d = 1.0 the former PDF changes into P(Jij) = \{${Jij — 
2) + 8(Jij)} implying that some of the bonds are missing, diluted bonds. If 
we also set c = 0.0 in Eq. (3), then some of the sites are occupied either by 
non magnetic particles or are empty, since the respective PDF (3) converts 
into P(hi) = p5(hi — hi) + q5(hi). These cases are shown in Fig 1(a) with each 
graph labelled by the respective site probability p with the graphs forming two 
groups: the upper one corresponding to c = d — 0.0 (site diluted system with 
vanishing exchange integral deviation) and the lower group corresponding to 
c = 0.0,d=1.0, site and bond diluted system. The choice d — 1.0 implies that 
the deviation is of the order of the exchange integral (d ~ \J\) so that some 
bonds may vanish according to the Eq. (2) whereas the remaining ones are 
strengthened. Also, other phase diagrams are presented in Figs. l(b,c) with 
each line labelled by a p- value. In Fig. 1(b), with d fixed as d = 0.5, the five 
upper solid lines correspond to c = 0.5, whereas the five lower dashed ones to 
c = 1.5. The topmost solid curve (i) is for p = 0.1 and as p increases towards 
p = 1.0, (curve (t>)) the respective curve lies below that corresponding to a 
lower value for p, the higher the p- value the lower the respective curve lies; this 



causes a gradual reduction of the ferromagnetic region. As far as the second 
group is concerned in Fig. 1(b) (dashed ones), an inversion of the order of the 
curves occurs in comparison to the upper solid ones with respect to p, that 
is, now the lowermost one (i) corresponds to p — 0.1, whereas as p increases 
towards to p = 1.0 (v ), the respective line lies above that with a lower value for 
p, so that, the higher the p- value the higher the curve lies, thus an enlargement 
of the ferromagnetic region takes place. However, before the inversion of the 
order of the lines takes place in panel (b), one observes coincidence of the 
phase-diagram boundaries for the all the p-values for d = 0.5, c = 1.0; in this 
case hi = \h 2 \. However, the order of the curves is unaltered if the temperature 
T is drawn with respect to the random field h\ for a fixed value of c (c = 0.5), 
whereas d varies as d = 0.00001 (solid curves) and d = 0.5 (dashed ones), 
Fig. 1(c), with the graphs for the same d— value forming a separate group, 
whereas in Fig. 1(b) they form a single group. A distinctive feature resulting 
from the plots in Fig. 1 is that all the lines of the phase diagrams coincide for 
small values of hi, beginning from the same point on the T-axis for hi = 0.0, 
implying that the critical temperature is the same regardless of c and p but 
depending only on d. Also, in this figure the critical temperature for hi =0.0 
is greater the smaller the d— value is. All the critical lines, separating the FM 
and PM phases, are of second order phase transitions. 
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Fig. 2. (Color online) The critical temperature T as a function of the exchange-in- 
tegral deviation d. In both panels, the upper solid graph (i) corresponds to hi = 0.0 
irrespectively of p (panel (a)) or c (panel (b). In panel (a), the lower graphs are for 
fixed c and hi as c = 0.1, hi = 1.0 and labelled by the selected p— values. In panel 
(b), the lower graphs are for fixed p and hi as p = 0.75, hi = 1.0 and labelled by 
the selected c— values. In both panels the higher the labeling value the smaller the 
extent of the ferromagnetic phase. 

The dependence of the critical temperature T c on the exchange-integral devi- 
ation d = A/ J appears in Fig. 2 for fixed random field hi = 1.0 as well as 
for the case hi = 0.0 for comparison. An overall feature of both panels is that 
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Table 1 



,(i) 



The critical field h\ c — values for constant values c = 0.5, c? = (second column) 



,(2) 



(:-!) 



and h\ c — values for c = d = 0.5 (third column) with respect the parameter p. h\ 
refers to the parameter d as a variable with constant values c = 0.5, p = 0.1 (fifth 
column), whereas h\ refers to the parameter c as a variable with constant values 
d = 0.5, p = 0.1 (seventh column). 



P 



h 



(i) 
l,c 



h 



(2) 
l,c 



(1 



h 



(3) 
l,c 



/?. 



(1) 
l,c 



0.0 6.736241 6.194181 0.0 5.904740 0.0 11.124309 

0.1 5.900871 5.406162 0.1 5.888921 0.1 9.678039 

0.2 4.421102 4.867123 0.2 5.840427 0.2 8.624253 

0.3 4.849948 4.435107 0.3 5.758102 0.3 7.543436 

0.4 4.488253 4.099344 0.4 5.640081 0.4 6.365760 

0.5 4.190953 3.816993 0.5 5.482039 0.5 5.482039 

0.6 3.940628 3.572388 0.6 5.277713 0.6 4.768053 

0.7 3.725971 3.363212 0.7 5.017263 0.7 4.189970 

0.8 3.538211 3.190681 0.8 4.686686 0.8 3.716970 

0.9 3.372639 3.053922 0.9 4.269706 0.9 3.324434 

1.0 3.224684 2.894779 1.0 3.768020 1.0 2.993567 

each line starts from a different point on the T— axis for d = 0.0 as compared 
to Figs. 1. In Fig. 2(a), c was fixed as c = 0.1 with the respective curves 
labelled by a p— value, whereas in Fig. 2(b), the constant quantity was p as 
p = 0.75 with the individual curves labelled by a c— value. In both panels the 
boundary line for hi = 0.0 lies above the ones for non zero values of hi, thus 
defining the widest ferromagnetic phase space in comparison to that for non 
zero h[s implying that the effect of including the random field is to reduce the 
ferromagnetic region; also the boundaries are arranged one below the other in 
ascending order either for p or c so that each time p or c is increased the re- 
spective ferromagnetic region is reduced in comparison to that with a smaller 
p or c. 

The critical field hi jC , the field for which the critical temperature vanishes, de- 
pends on the thermodynamic route followed, since the current random system 
contains various parameters, namely, c, d, p; consequently, in order to estimate 
hi fi , two of these parameters are fixed, only the remaining one varies. Initially, 
c and d are kept fixed as c = 0.5, d = and later c — d — 0.5; the values of the 
respective critical field hi iC appear in Table 1 for selected p— values for both 
choices. In the former case, these values are fitted by the sixth-degree poly- 
nomial P(x) = 6.73612 - 9.97119a; + 19.00728a; 2 - 30.76579a; 3 + 33.56416a; 4 - 
20.72121a; 5 + 5.37542a; 6 , whereas in the latter case by the fifth-degree polyno- 
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mial P(x) = 6.19323 - 9.39129a; + 18.25539a; 2 - 27.81661a; 3 + 23.72744a; 4 - 
8.07256a; 5 , with x = p in both cases. Also, in the same Table, there appear the 
values for the critical field h[ c (fifth column) corresponding to the constant 
quantities c = 0.5, p = 0.1 with parameter d as a variable, for which the best 
fit polynomial is P(x) = 5.90487 - 1.6687a; 2 + 0.39594a; 3 - 0.8657a; 4 (x = d); 
moreover, in the seventh column there appear the values for the critical field 
h\ c corresponding to the constants d = 0.5, p = 0.1 and parameter c as a vari- 
able, with best fit polynomial P(x) = 11.12544 - 14.35749a; + 6.410899a; 2 - 
0.155474a; 4 (x = c). 



M 




Fig. 3. (Color online) The thermal dependence of magnetization M(T,hi), order 
parameter. In panel (a), the fixed parameters are d = 0.1, p = 0.1, hi = 1.0; each 
graph is labelled by the random field ratio c = hilh\\ c = 0.1 graph (i), c = 1.0 graph 
(ii) and c = 2.0 graph (iii). In panel (b), the fixed parameters are c = 3.0, d = 0.1; 
each graph is labelled by p and hi as: innermost graph p = 0.1, hi = 0.2, middle 
one p = 0.6, hi = 0.2 and outermost p = 0.1, hi = 0.0. 

In addition to the phase diagram, another important quantity is the mag- 
netization; its thermal behavior appears in Fig. 3 for various values of the 
system's parameters c,d,p,hi, resulting by solving numerically the Eq. (20). 
In panel (a) the parameters d,p,hi are fixed as d — 0.1, p = 0.1, hi = 1.0; 
the individual graphs are labelled by the random field ratio c = ^2/^-1, the 
outermost graph (i) corresponds to c = 0.1, the middle one (ii) to c = 1.0 
and the innermost one (iii) to c = 2.0. The magnetization exhibits its normal 
behavior but it is stronger for the smaller value of this parameter c = 0.1 
(for c = 0.0 the respective graph is very close to that for c = 0.1 due to 
the finiteness of the calculations), which is in accordance with the respec- 
tive plot in Ref. [12] although in the latter publication the random field is 
not included. In the latter case (c = 2.0) the stronger negative random field 
/12 = chi yields smaller magnetization because this field competes strongly the 
first term in the Hamiltonian (1), which favors the parallel orientation of the 
spins. In panel (b) the parameters c, d, hi are fixed as c = 3.0, d — 0.1, hi — 0.2, 
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the plots are labelled by the site probability p = 0.1 (innermost graph) and 
p = 0.6 for the middle graph; for comparison the case for the lack of the ran- 
dom field was also included, which is the outermost graph corresponding to 
c = 3.0, d = 0.1, p = 0.1, hi = 0.0. For small values of the temperature, the 
magnetization values are constant, independent of the relative parameters. 

A susceptibility-like quantity Xt(^-i) can also be calculated as the first order 
derivative of the magnetization with respect to the applied external random 
magnetic field, namely 



Xr(hi 



dm dm dh 1 dm 
dhi dho dh\ J dho 



(23) 



so that 



JXrih) 



dm 4mf^ + 4m 3 f^ 
dh 1 - 4Qi4 - 12m 2 Q 34 



(24) 



or, for calculational purposes, 

/ T ,, s\-* 1 - 4g M - 12m 2 Q 3 4 



4m^g^ + 4m 3 

oho 



9Q34 

dho 



(25) 



since h = J hi, although the product (J hi) was earlier set as hi, now we 
divert for calculational purposes. 




Fig. 4. (Color online) The inverse magnetic susceptibility-like vs temperature T, 
with constant c = p = hi = 0.1; graph (i) corresponds to d = 0.1 and graph (ii) to 
d= 1.0. 

The temperature dependence of the inverse susceptibility-like is shown in 
Fig. 4 for selected values of the parameters c, d,p, hi. In both graphs, the sus- 
ceptibility diverges twice, once at the zero temperature (T = 0, finite clusters 
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contribution) and once again at the respective critical point T c , infinite clus- 
ter contribution. Although both graphs increase smoothly for small T— values, 
they decrease steeply for larger ones. A similar behavior for the susceptibility 
was also found by Kaneyoshi et al [24]. 

4 Conclusions and discussions 

In the current investigation we have determined the phase diagrams, critical 
temperatures, magnetization profiles and susceptibility of the nearest-neighbor 
spin-1/2 Ising model on a square lattice, in the presence of competing random 
fields as well as fluctuating random bonds via EFT, based on the differen- 
tial operator technique and Callen identity. The topology of the system is 
taken into account through the coordination number. The EFT method is an 
improvement on the MFA, since it provides vanishing critical temperature for 
one-dimensional systems, also for two-dimensional systems the respective crit- 
ical temperature is closer to the exact one than that of the MFA. A variety of 
phase diagrams were obtained with respect to the various system's parame- 
ters, namely, c, d, hi, p. The extent of the ferromagnetic region does not show 
systematic variation on varying the system's parameters; it varies non mono- 
tonically for constant d and varying c, whereas in other circumstances it varies 
monotonically. The magnetization profile was determined with respect to tem- 
perature as well as the random field hi. An immediate result of the presence 
of the random field is to reduce the extent of the ferromagnetic region, in gen- 
eral, as well as numerical value of the magnetization as c increases because of 
the strengthening of the random field h 2 = chi and the subsequent increased 
tendency of this field to prevent spins to attain the parallel configuration and 
the competition with the other random field h\. 

Although the EFT technique improves on the MFA, it presents the same short- 
comings in the critical region (non-classical region) concerning the critical ex- 
ponents like the MFA and Landau theory, because the fluctuations occurring 
in the critical region, as the transition temperature is approached, become im- 
portant and the non-classical behavior is observed; these fluctuations are not 
taken under consideration properly by EFT or MFA, thus a relative criterion, 
called Ginzburg criterion, determines how close to the transition temperature 
the true critical behavior is revealed [25]. This criterion relies on any thermo- 
dynamic quantity but the specific heat is usually considered for determining 
the critical region around T c where the mean field solution cannot describe 
correctly the phase transition. The MFA is valid for lattice dimensionality 
greater than the upper critical dimension d u = 4 in case of presence of only 
thermal fluctuations. However, in the current case the presence of random 
fields enhances fluctuations causing the critical region to be wider than the 
one due only to the thermal fluctuations and the upper critical dimension is 
increased by 2 to d u = 6 [26,27]. Occasionally the non classical region is ex- 
tremely narrow so that the respective critical behavior expected by Landau 
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or MFA is observed because the fluctuation region is very narrow and hardly 
accessible for experimental observation; such a system is the weak-coupling su- 
perconductor in three dimensions for which the respective non classical region 
is l^cfll < 1CT 16 (tcR is the reduced temperature). However, on reducing the 
space dimension as in the case of the weak-coupling superconductor in two di- 
mensions, the critical exponents have their classical values up to \tcR.\ = 10 -5 , 
thus the reduction of the space dimensionality causes serious repercussions on 
the behavior of the physical system; on the contrary, for the superfluid helium 
transition the classical region extends up to |£or| < 1-0 so that fluctuations 
are detectable [28,29,30,31]. In addition to superconductivity, the extent of 
the non classical region for the ferroelectric system triglycine sulfate (TGS) is 
relatively small and its critical exponents have the respective classical values 
up to \t CR \ = 1.5 x 1(T 5 [32,33,34]. 

The results obtained in the current investigation by using the EFT provide a 
basis for a comprehensive analysis by more sophisticated methods. However, 
they are of no less importance, since they show, nevertheless, the expected 
phenomena to be observed. 
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